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^ Abstract 

!• We calculate the magnetic quadrupole moment (MQM) of the deuteron at lead- 

• ^H ing order in the systematic expansion provided by chiral effective field theory. We 

/\ take into account parity {P) and time-reversal (T) violation which, at the quark- 

^ gluon level, results from the QCD vacuum angle and dimension-six operators that 

originate from physics beyond the Standard Model. We show that the deuteron 

MQM can be expressed in terms of five low-energy constants that appear in the 

P- and T-violating nuclear potential and electromagnetic current, four of which 

also contribute to the electric dipole moments of light nuclei. We conclude that 

the deuteron MQM has an enhanced sensitivity to the QCD vacuum angle and 

that its measurement would be complementary to the proposed measurements of 

light-nuclear ED Ms. 



Permanent electric dipole moments (EDMs) of particles, nuclei, atoms, and molecules 
are powerful probes for physics "beyond" the Standard Model (SM) of particle physics. 
EDMs violate both parity (P) and time-reversal (T) invariance (or, if we invoke the CPT 
theorem of relativistic quantum field theory, CP invariance). EDMs could be low-energy 
manifestations of some source of P and T violation {f^) that originates at an energy scale 
comparable to or even higher than that accessed by the LHC At the current experimental 
accuracy, electroweak quark mixing can be ignored, and the only Standard Model (SM) 
source that can impact EDMs is the QCD vacuum angle, i9 [1]. A nonzero measurement 
of a hadronic EDM could be due to either beyond-the-SM f^ sources or a finite value of 
9, even though the current upper bound on the neutron EDM already limits this angle 
to a very small value. Among the possible beyond-the-SM f^ sources, those represented 
by operators of effective dimension six [2j are expected to dominate: the quark EDM 
(qEDM), the quark and gluon chromo-electric dipole moments (qCEDM and gCEDM, 
respectively), and two independent four-quark (FQ) interactions. From a theoretical 
point of view, it is clearly an important priority to disentangle the dimension-four SM 
^-term and the non-SM sources [3]. 

The deuteron is an excellent candidate for a sensitive EDM search in a storage ring |1]. 
In a recent paper [5] we addressed the EDMs of light nuclei, including the deuteron, within 
chiral effective field theory (EFT). The major advantage of such an EFT approach is its 
direct link to QCD in that it exploits the different chiral properties of the fundamental 
f^ sources. Moreover, the power-counting scheme allows for a controlled framework 
such that the theoretical uncertainties can be estimated and the results can be improved 
systematically. We showed that the EDMs of light nuclei can be expressed in terms 
of essentially six f^ parameters, or low-energy constants (LECs). These LECs can in 
principle be calculated from the underlying f^ sources by solving QCD at low energies, 
in particular by lattice simulations. Lacking that, the size of the LECs can be estimated 
by naive dimensional analysis. We concluded that the EDMs of various light nuclei can 
give crucial complementary information about the fundamental f^ sources. 

Since it is a spin-1 particle, the deuteron has one other static f^ electromagnetic mo- 
ment, the magnetic quadrupole moment (MQM). In this paper, we address the deuteron 
MQM in the same framework that was used in Ref. |3]. It was shown in Ref. [31| that, in 
addition to LECs that contribute to the EDMs of light nuclei, the deuteron MQM depends 
also on f^ pion-nucleon-photon interactions. Moreover, it was argued that only for the 9 
term is the deuteron MQM expected to be larger than the deuteron EDM (in appropriate 
units). For the beyond-the-SM sources, the MQM is expected to be of similar size or 
somewhat smaller than the EDM. This indicates that a measurement of the deuteron 
MQM, if possible, could play a central role in separating the various f^ sources. 

The conclusions of Ref. [3] were based on a chiral EFT in which the one-pion exchange 
nucleon-nucleon (AW) force is treated in perturbation theory [6]. This approach, expected 
to be valid at low energies, allows one to give analytical results for the deuteron EDM 
and MQM, which, moreover, can be extended to sub-leading orders. On the other hand, 
in NN scattering it was found that the results do not converge in certain partial waves 
for momenta somewhat lower than expected [7j. It is therefore important to check the 



results of Ref. [3] in the framework of a chiral EFT that treats pions nonperturbatively [8]. 
For the deuteron EDM it was found that the two EFTs gave similar results [5]. In this 
article we investigate this for the deuteron MQM as well and we address the question, 
to what extent a possible measurement of the MQM could be of help to separate the 
different fj" sources. We also compare our results to previous studies of the deuteron 
MQM [HI Hn] • In particular, Ref. [TU] used traditional meson-exchange NN models and a 
general f^ NN interaction [TOl [H]. We use the codes of Ref. [I0|, but adapt and extend 
the framework (the f^ NN potentials and currents) to chiral EFT with nonperturbative 
pions. We follow the hybrid approach of Ref. [5] in which nonperturbative pion exchange 
is embedded into the P- and T-conserving {PT) wave functions of modern, "realistic" 
potentials P^ [T3] . In principle our framework can be applied to the calculation of the 
MQM of other light nuclei as well. 

We start our discussion of the deuteron MQM from the effective hadronic Lagrangian 
involving the low-energy degrees of freedom: nucleon {N), photon (A^), and pion (tt). The 
PT part, which originates from the quark kinetic (color- and electromagnetically gauged) 
and mass terms in QCD, is well-known [TU [S]. At leading order (LO) it consists of the 
standard pion- nucleon axial- vector coupling, qa = 1.27, and the pion-nucleon-photon in- 
teraction obtained from gauging the qa term. The pion-photon interactions stem from 
the pion charge. At next-to-leading order (NLO), the photon couples to the nucleon via 
the covariant derivative in the nucleon kinetic term and via the isoscalar and isovector 
magnetic moments, respectively Kq = —0.12 and Ki = 3.7. The f^ effective Lagrangian 
results from the most general /^ Lagrangian at the quark-gluon level. It was derived 
and discussed in detail in Refs. [I5l HE]. The key observation is that the different fun- 
damental fl" sources transform differently under chiral symmetry and therefore generate 
different hadronic interactions. Given enough independent observables it will be possible 
to disentangle them [3]. 

We present here only the terms that are relevant for the LO MQM calculation, which 
depends on five fl" interactions. The Lagrangian 

Crr = -^N{goT-7r + g^7C3)N+^e^''''^Vo.NS^T-7rNF,, 

-^ TV -^77 

+CiNNd^{NS^N) + C2NTN ■ V^{NS^tN) , (1) 

contains isoscalar (^0) and isovector {cji) nonderivative pion-nucleon couplings, an isoscalar 
pion-nucleon-photon coupling (ct^), and two short-range f^ NN interactions {Ci, C2)- 
Here F^^ = 185 MeV is the pion decay constant, r are the Pauli matrices in isospin space, 
Vfi (Sfj^) is the nucleon velocity (spin), e'^^°''^ (with e°^^^ = 1) the completely antisymmetric 
tensor, F^,^ = d^A^ - dyA^ the photon field strength, and {V^)ab = Sabdf, + ee^^abA^ + 
0{-K ■ df^Tz/F^) a chiral covariant derivative. Except for c^, these interactions play a role 
in the calculation of nuclear EDMs [21 El HTj as well. The scaling of the LECs in terms of 
the pion mass (m^r), the characteristic QCD scale (Mqcd ~ '^'^P-k), and the scale oi f^ 
physics beyond the SM is given in Ref. ^ for go, gi, Ci, and 6*2. The pion-nucleon-photon 
coupling Ctt has the same scaling as the short-range contribution to the isoscalar nucleon 
EDM, do t3j. 



The calculation of the deuteron MQM can be divided into two contributions. The first 
contribution conies from an insertion of the f^ electromagnetic two-body current Jfj^. 
The current has to be two-body since the constituent nucleons, being spin-1/2 particles, 
do not possess a MQM. The second contribution comes from the electromagnetic current 
JpT upon perturbing the wave function of the nucleus with the f^ potential \^, such 
that the wave function obtains a f^ component. This current can be one- or two-body. 
The required f^ potential Vfrp and current Jfrp can be calculated from Eq. (Il| . 

To first order in the f^ sources, the deuteron MQM is thus a sum of two reduced 
matrix elements, 

Md = ^{{^d\\M\\m,)+2{^!a\\M\\^S) . (2) 
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The deuteron ground state \^d) and its parity admixture \^ d) are the solutions of homo- 
geneous and inhomogeneous Schrodinger equations, 

{E-Hpt)\^j) = 0, (3) 

{E-Hpt)\^,) = %|*,), (4) 

respectively, where Hpx is the PT Hamiltonian. The MQM operators M and M are 
obtained from the corresponding currents JpT and J^, respectively. The Cartesian com- 
ponent along the z direction, M33, which is proportional to the spherical harmonic ¥2, 
takes the form 

M33 = 2 / Sx X3 [x X J(f))3 , (5) 



where x is the position where the current density is probed. Given the current in momen- 
tum space, 

/(a^/.'.e-<-V>-). (6) 

M33 can also be derived as 

M33 = -2 hm (V,3V,, J2(g1 - V,3V,, Ji(g-)) . (7) 

i3-s>0 

The three classes of contributions to the deuteron MQM described above are shown 
in Fig.[T} In order to decide which diagrams give the main contribution to the MQM we 
apply the power counting rules outlined in Ref. |5]. These rules provide an expansion in 
Q/Mqcd, where Q is the generic momentum in the process, in this case of the order of the 
deuteron binding momentum. We use Q ~ m^ ~ Fj^, as standard in xPT. For each class 
we take the PT vertices from the interactions described above and the f^ LO vertices 
from Eq. ([I|. The iteration of the LO PT potential is not suppressed, and is necessary 
among nucleons in reducible intermediate states, as indicated in diagrams (b) and (c) of 
Fig.fl] Such iteration among nucleons before and after any f^ insertion builds up the 
PT wave function, represented in Fig.[T]by the triangles. In the following power counting 
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Figure 1: The three general classes of diagrams contributing to the deuteron MQM de- 
scribed in the text. Solid and wavy lines represent nucleons and photons. The large 
triangle denotes the PT wave function and the oval iterations of the PT potential. The 
nucleon-photon vertex is the nucleon magnetic moment or the convection current, and 
the photon attached to the oval is the PT two-body current. The black squares without 
(with) an attached photon represent the f^ potential (two-body current). 



we omit this overall factor. The scaling in terms of the LECs in Eq. ([I]) of diagram (a) is 
then 



Da = e 



F^MqcbJ ^^ \ F^ Mqcd 
while diagrams (b) and (c) scale as 
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QCD 



(9) 



Which contribution dominates depends on the fundamental fj" source, the dimension- 
four ^-term and the dimension-six terms: qEDM, qCEDM, gCEDM, and FQ. 

• For the 6 term, only the isoscalar coupling qq plays a role at LO. In order to generate 
gi, the 6 term requires an insertion of the quark mass difference, which causes a 
relative suppression of gi relative to ^o by a factor em^/Mqf^j-, [TJ], where e = 
{rrid — mu)/{rnu + rrid)- The other couplings, c^r and Ci^2, also contribute at sub- 
leading orders. Since ^o enters in principle through the three classes of diagrams 
with similar factors, all these classes are equally important. 

• For the qCEDM, we need both pion-nucleon interactions, since there is no relative 
suppression of gi. Again, all three classes are, a priori, equally important. 

• For the qEDM, the purely hadronic interactions are suppressed by factors of the 
fine-structure constant, and only c^ is important. Thus only diagram (a) matters. 

• For the chiral-invariant (%!) f^ sources (gCEDM and FQ), the pion-nucleon interac- 
tions, which break chiral symmetry, are suppressed by a factor m^/Mq^j-, compared 
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Figure 2: Diagrams contributing to the/7^ two-nucleon current. Solid, dashed, and wavy 
lines represent nucleons, pions, and photons. A square marks a f^ interaction and the 
other vertices PT interactions. Only one topology per diagram is shown. 



to the f^ short-range AW interactions, which conserve chiral symmetry. The latter 
are therefore as important as the pion-nucleon interactions. For the xl sources again 
all three diagrams could be important. 

We now turn to the calculation of the various ingredients in these diagrams, starting 
with the currents. For the nucleons we use incoming momenta pi = P/2 + p and p2 = 
P/2 — p, and outgoing momenta Pi = P /2 + p' and P2 = P /2 — p'. The photon 
momentum q = P — P is outgoing. For convenience we introduce k = p — p'. The spin 
(isospin) of nucleon i is denoted by a'^^V^ ('''^^V^)- 

The relevant f^ two-body currents that appear in diagram class (a) of Fig.jl] are 
shown in Fig. 12} Since the deuteron wave function is isoscalar, the f^ currents need to 
be isoscalar as well in order to contribute to the MQM. Only the current in diagram (c), 
which stems from qEDM, meets this requirement, and we find 
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. (10) 



In diagram classes (b) and (c) of Fig.fllthe photon interacts instead with a PT current. 
The PT one-body current in diagram (b) is either the nucleon magnetic moment or the 
convection current coming from the nucleon kinetic energy. 
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where i is the index of the nucleon that interacts with the photon. In diagram (c) we 
require the PT two-body currents depicted in Fig.|3j 
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With the methods outlined in Refs. p^ |5], these currents can straightforwardly be 
Fourier transformed to coordinate space, where we denote by a;'^*-' the position of nucleon 
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Figure 3: Diagrams contributing to the PT two-nucleon current. 



and by r = x^^^ — x^"^^ the relative position. We also introduce the Yukawa function 
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The MQM operators are, from Eq. (10), 

= 2^T«-T(2)[(3a«4^)-a«-x(i 



M 



)W(2) 



(2)^(2) _ -,(2) . -j(2)A -.(1) 



3 (To X 



3 -^3 



a^ ' ■ x^ ' ] a 



VrU{r) 



(14) 



and from Eqs. (11) and (12) 
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where p^*^ = — zV^-ci). 

The last ingredients we need are the deuteron wave function ^ ^ and its parity admix- 
ture ^d- They have been calculated in Ref. [TU] using modern high-quality phenomeno- 
logical PT potentials and a f^ potential dominated by pion exchange and extended with 
heavy-meson exchange. In the EFT spirit, it would be best to calculate '^d from a PT 
Hamiltonian fully consistent with chiral symmetry and renormalization-group invariance, 
but an accurate fit to two-nucleon data with these properties does not yet exist. For 
those sources for which the MQM is dominated by long-distance physics, we expect no 
significant differences when using the phenomenological potentials, and thus for \E'^ we 
use the results of Ref. [10]. For comparison with Ref. [5], we give below numbers cor- 
responding to the Argonne wis (AV18) potential [12]. Differences with results obtained 
from the Nijmll and Reid93 potentials [13] are within a few percent, except for the cases 
of C*! 2 insertion — for which more details will be provided later. 

The ground state of the deuteron is mainly a ^Si state with some ^Di admixture. The 
matrix element of M is found to be 



^ (^rf||M||^rf) = 0.07— efml 



(16) 



The relevant matrix elements of M are obtained with "^d from the LO f^ two-nucleon 
potential. The general fl" NN potential was derived in Ref. p3] and we summarize the 
relevant parts here. In coordinate space the potential is given by 
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(17) 



In this expression, at LO go originates from 6 term, qCEDM, and %! sources; cji from 
qCEDM and xl sources; and Ci from xl sources only. 

If the parity admixture comes from a gi pion exchange, the deuteron wave function 
acquires a ^Pi component. In order to get back to the ground state, the deuteron can 



couple to the photon via the isovector one-body or the two-body currents in Eqs. (|ll|) 

r 



and (12), respectively. The result is 
2 
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0.031(1 + /€i) + 0.003 + 0.008 



91 



efm 



where the first and second terms come from the isovector magnetic moment due to the 
spin and convection current, respectively; and the third term from the PT two-body 
current. 

After a go pion exchange or an insertion of Ci^2 the deuteron wave function obtains a 
^Pi component instead. In this case, the deuteron needs to couple to the isoscalar nucleon 
magnetic moment or an isoscalar two-body current, which is not present at LO. The result 
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For the contact interaction with 



Co = Ci — 3C2 



;i + fi:o)efm^ (19) 



(20) 



we apply a strategy followed in Refs. [20l |5]: it is simulated in our calculations by a 
fictitious heavy-meson (of mass m) exchange, since 



m^Co 
Airr 
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as m goes to infinity. As shown in Fig.|4| when m reaches 2.5 GeV, the results converge 
at about < 10% level, so we report the above numbers at this scale. While in this figure 
one sees good consistency between the AV18 and Nijmll results, the Reid93 result is off 
by a factor of 2. The main reason is that the Reid93 potential generates a deuteron S 
state whose short-distance wave function is enhanced, leading to more sensitivity to Co- 
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Figure 4: Deuteron MQM due to the F^(— C1+3C2) short-range /^ interaction as function 
of a regulating mass m, for various PT Hamiltonians. 



The large discrepancy between different potentials suggests that for xl sources, for which 
Cq contributes to the MQM at leading order, a fully consistent calculation of \E'd within 
EFT is necessary if this part of the matrix element needs to be known better than within 
a factor of 2. 

The dependence of the deuteron MQM on (70,15 Co, and Ct, was studied in Ref. [3] in 
a framework where pion exchange is treated perturbatively. At LO in that framework, 
the coefficients in front of go{l + kq)/Ft, and gi{l + Ki)/ F^, were found to be, respectively, 
—0.146 and —0.049, in agreement with the results in Ref. [9j where a zero-range approx- 
imation for the NN interaction was assumed. Considering the large intrinsic uncertainty 
(~ 30%) in the perturbative-pion calculation, the perturb at ive-pion gi coefficient is in 
reasonable agreement with Eq. (18). A similar agreement was found for the deuteron 
EDM [3]. On the other hand, the perturbative-pion g^ coefficient is three times larger 



than Eq. (19), suggesting that the effects of additional pion exchanges, neglected in the 
LO perturbative calculation, are larger in the ^Pi channel than in the ^Pi channel. We 
have verified that, if the tensor force is ignored and the same strong force is assumed for 
both the ^Pi and ^Pi channels, the ratio of the strong parts of the MQM matrix elements 
due to the isoscalar and isovector one-body currents becomes 3, which is consistent with 



Refs. [Oils]. Preliminary results of an NLO calculation in the perturbative-pion framework 
indicate that, indeed, NLO corrections influence the go coefficient by a larger amount than 
the gi coefficient [2T] . 

In the framework of perturbative pions, the PT convection and two-body currents in 
diagram (b) and (c) of Fig. [I] enter at NLO and are expected to be smaller than the 
contribution from the isovector magnetic moment. This agrees with the numerical results 



in Eq. (18) where the convection and two-body currents only enter at the ~ 5% level. 
These currents would have been small (~ 30%) even if the isovector magnetic moment 
had been more natural, that is, if 1 + /ti were ~ 1. 



The result for c^ in Eq. (16) is somewhat smaller than expected from the power 
counting estimate in Eq. (IS]), A^d(qEDM) ~ 0.2F^C7rfm^, and is more in line with 
the expectation of the perturbative-pion calculation, 0{'yc.n-/MNNMQCD) ~ 0.07F.„CTrim.^, 
where M^n = 4:7iF^/g\mN is the characteristic scale where pions become nonperturba- 
tive [6] and 7 the deuteron binding momentum. A more detailed comparison with the 
perturbative-pion calculation is complicated by the appearance of a short-range current, 
which is needed for renormalization purposes. Neglecting the counterterm and using for 



the renormalization scale fx = Mnn, the perturbative result becomes O-OSF-^CT^iT^ , in 



good agreement with Eq. (16). A comparison between the contributions from the short- 
range f^ NN interactions in the perturbative and nonperturbative calculations is not 
useful because the LEG Cq includes different physics and thus has different scalings in the 
two EFTs. 

The deuteron MQM was previously calculated in Refs. P ^0\, in which the deuteron 
MQM was assumed to be dominated hj fl" one-pion exchange (OPE). Since these calcula- 
tions did not use the chiral properties of the fundamental f^ sources, the /^ pion-nucleon 
interactions were assumed to be all of the same size. Our analysis shows that these as- 
sumptions only hold in case of a qCEDM. For g^ and gi OPE we confirm the results in Ref. 
[TU] . but our go result is 3 times smaller than the result in Ref. [S] due the discrepancy 
discussed above. 

We are now in the position to discuss the results for the various f^ sources. It was 
noted in Ref. p] that the observation of a deuteron MQM, together with the nucleon 
and deuteron EDMs, could provide important clues to separate the various f^ sources. 
In particular, it was concluded that only for chiral-symmetry-breaking, but isoscalar, 
sources like the QCD 6 term is the deuteron MQM, in appropriate units, substantially 
larger than the deuteron EDM. For chiral- and isospin-breaking sources, like the qCEDM, 
the deuteron MQM and EDM are expected to be of the same size, and a measurement 
of both would fix the couplings go and gi, allowing a prediction of other f^ observables, 
like the ^He EDM. For the %! sources and the qEDM, in the perturbative-pion approach 
the MQM depends on one- and two-body LECs that do not contribute to the EDM. 
For these sources the MQM was found to be of the same size, or slightly smaller, than 
the EDM, but an observation of the MQM would not give us predictive power. These 
conclusions, based on the perturbative-pion power counting, are confirmed here by the 
nonperturbative results. 

For the QCD 6 term, the deuteron EDM is dominated by the isoscalar nucleon EDM, 



dd ~ 2(io [3j ^. A naturalness lower bound on the isoscalar nucleon EDM is provided by 
the non-analytic terms stemming from the pion cloud [22j, \do\ ^ .01{\go\ / F^^) efm. The 
deuteron MQM is dominated by the go piece in Eq. (19). Combining the two with the 
deuteron mass m^, we find 
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(22) 



At LO in the perturbative-pion approach, this ratio is about three times larger, as dis- 
cussed above. Nonetheless, the nonperturbative calculation confirms that for isoscalar 
chiral-breaking sources the deuteron MQM is expected to be larger than the EDM in 
units of rrid. 

In case of the qCEDM, where ^o and gi have similar scalings, both the deuteron EDM 
and MQM are dominated by pion exchange. At LO, the EDM depends on cji only |5], 
and the MQM on the ^o,i contributions in Eqs. (|I8|) and (|l9|). The MQM/EDM ratio 
becomes 



"^ ^ 1.6(1 + «:i) + 2.2(1 + «:o)^ + 0.6, 

dd gi 



(23) 



which formally is 0{1). However, due to the large anomalous isovector magnetic mo- 
ment, numerically the ratio could be (9(10). This means that the measurement of a 
large MQM/EDM ratio does not necessarily imply that the 6 term is the dominant f^ 
mechanism. 



If a qEDM is the dominant f^ source, the MQM is given by Eq. (16). It is solely 
coming from a two-body f^ current. The EDM is given by the sum of the neutron and 
the proton EDM, dd — 2do [3]. Combining these results gives 



mdMc 
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0.7- 



4' 



(24) 



which is 0{1) by naive dimensional analysis (NDA). As we observed in the previous 
discussion, the matrix element of the operator with coefficient c.^ is smaller than the 
power counting estimate, so a more accurate conclusion is that, for f^ from the qEDM, 
the deuteron MQM is expected to be slightly smaller than the EDM, in agreement with 
Ref. [3]. 



Finally, for chiral-invariant sources, the MQM is dominated by the sum of Eqs. (18) 
and (|19|). From the NDA estimates in Ref. [5], we find F^Co/go = 0{F^/ml) ~ 2, 
such that the Cq contribution only enters at the ~ 10% level. The deuteron EDM at 
LO formally depends on gi and the isoscalar nucleon EDM, but numerically the latter is 
expected to dominate, with pion-exchange corrections at the ~ 15% level [5]. Ignoring 
the numerically small convection and two-body currents. 



mdMc 
dd 



0.21(1 + /to) 



9q_ 

fA 



0.15(1 + /ti) 



7r"0 



9l_ 

Ft, do 



efm. 



(25) 
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By power counting we would expect the ratio to be 0{1), but since the deuteron is 
weakly bound, pion exchange is smaller than expected. Using the NDA estimate |(io| ~ 
5(|(7o,i|/-^7r) efm, we conclude that for %! sources the deuteron MQM should be smaller 
than the EDM. 

In conclusion, we computed the deuteron MQM for various fl" sources: the QCD 6 
term, quark EDM, quark and gluon chromo-EDMs, and chiral-invariant four-quark op- 
erators. We performed these computations at leading order in the framework of chiral 
EFT, with pions treated nonperturbatively. The same parameters as in the correspond- 
ing calculation of light-nuclear EDMs appeared here, except for the quark EDM, which 
involves an independent short-range two-nucleon current. While the results confirm the 
qualitative conclusions of Ref. [3J, there are important quantitative differences. Due to its 
enhanced sensitivity to the QCD 6 term, a potential measurement of the deuteron MQM 
would be complementary to one of the deuteron EDM. 
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